D*Dp vertex from QCD sum rules 
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, constant: go^op = 4.1 ± 0.1 GeV~^. 
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I. INTRODUCTION 



Over the last years the strong interaction of charmed hadrons among themselves and with other species of hadrons 
rS I has received increasing attention. From the discovery of charmed mesons in the seventies until the late eighties there 
' was no motivation to study in detail the interactions of these particles. In the early nineties there was a series of 
papers trying to compute the cross section of a J/ip with ordinary light hadrons. The motivation came from the heavy 
(— I ' ion program running at CERN and later at RHIC. At that time J/ip suppression was considered as a signature of 
I— ' quark gluon plasma (QGP) formation and it was very important to know as accurately as possible the purely hadronic 
^ (non-QGP induced) charmonium suppression, which would be the background for the QGP signal. From 2000 on, it 
^ , became slowly clear that the physics of J/ip is much more complex than thought before and its simple suppression 
■ was no longer considered as a QGP signal and the subject lost interest. On the other hand, during those years, at the 
' B factories the collaborations BABAR and BELLE started to produce results. One of the important decay channels 
of the B mesons is into J/ip (plus other things). Moreover these collaborations found new charmonium states (the 
C^l X, the F's and the Z), which also decay into J/ip or into tp' . It has been conjectured that both B and the new 
charmonium states very often decay into an intermediate two body state with D^s and/or £'*'s, which then undergoes 
final state interactions, with the exchange of one or more virtual mesons. In order to calculate the amplitudes of these 
. processes we need to know the relevant vertices involving the charmed mesons. As an example of specific situation 
• where a precise knowledge of the D*Dp form factor is required, we may consider the decay X(3872) J/tp + p. 
^ \ As suggested in [l|, this decay proceeds in two steps. First the X decays into a. D - D* intermediate state and then 
JJ"! ■ these two particles exchange a D* producing the final J/ip and p. This is shown in Fig. lb and If of [ij. In order 
rN [ to compute the effect of these interactions in the final decay rate we need the D* Dp form factor. More generally, we 
^ i need to know all the charm form factors to correclty calculate the interaction of J/ with light hadrons and the final 
■ " " ' state interactions in B decays. These form factors have been calculated in the framework of QCD sum rules (QCDSR) 
[3| techniques in a series of works on vertices involving charmed mesons, namely D*Dn SH, DDp\^, DDJ/^p [6|], 
D*DJ/^ !7:], D*D*-K [1,11, D*D*J/iP [3, D,D*K, D^DK [HI , DDuj [13 and D* D* p^. 

In the present paper we calculate the D*Dp form factor with QCDSR. In the next section, for completeness we 
describe the QCDSR technique and in section III we present the results and compare them with results obtained in 
other works. 

II. THE SUM RULE FOR THE D*Dp VERTEX 

Following our previous works and especially Ref. [3|, we write the three-point function associated with the D*Dp 
vertex, which is given by 
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for an off-shell D meson, and: 

^^^}{P,P') - / d^^d^y e^P'-e^(p'-rt-^ (0|r{j^(x)j^(y)j-,^*^ (0)}|0) (2) 

for an off-shell p meson. The general expression for the vertices ([T]) and ^ has only one Lorentz structure. Equa- 
tions (tlj and ([2]) can be calculated in two diferent ways: using quark degrees of freedom -the theoretical or QCD 
side- or using hadronic degrees of freedom -the phenomenological side. In the QCD side the correlators are eval- 
uated using the Wilson operator product expansion (OPE). The OPE incorporates the effects of the QCD vacuum 
through an infinite series of condensates of increasing dimension. On the other hand, the representation in terms of 
hadronic degrees of freedom is responsible for the introduction of the form factors, decay constants and masses. Both 
representations are matched invoking the quark-hadron global duality. 



A. The phenomenological side 

The D*Dp vertex can be studied with hadronic degress of freedom. The corresponding three-point functions, 
Eqs. ([T]) and ([2]), are written in terms of hadron masses, decay constants and form factors. This is the so called 
phenomenological side of the sum rule and it is based on the interactions at the hadronic level, which are described 
here by the following effective Lagrangian [T5| 

Cd'Dp = - go- Dp e^*"^ ( D d^ps d^D; + hc^ (3) 

from where one can extract the matrix element associated with the D* Dp vertex. In the above expression we have 
g0i23 = Saturating Eqs. ([T]) and Q with the appropriate £>, D* and p states and making all the contractions we 
arrive at: 

rW(p,p') ^ A(fi(p^y^g2)e„^,.pV'' +h. r. (4) 
where h. r. means higher resonances and q = p — p' . The invariant amplitude A^*^) is given by 



for an off-shell D meson and 



\Hen - 9D>Dp{Q ) ^p, ^ ^ ^2)(p,2 ^ 

for an off-shell p meson. In the above expressions ~ —p^, P'^ = —p'^ and — —q^- 

The meson decay constants appearing in the equations above are defined by the vacuum to meson transition 
amplitudes: 

(Olj^ll?) = ^ , (7) 
rric 

and 

{V{p, e)\ji\0)^ my fve*^, (8) 

for the vector mesons V ~ D* and V = p. The form factor which we want to estimate is defined through the vertex 
function for an off-shell p meson: 

(i^*(p,A)|D(p')p(g,A')) =*gl3'^,,p(g')e"''^'e^(p)e^'(gH9. ■ (9) 

where (p) and (q) are the polarization vectors associated with the D* and p respectively. An analogous expression 
holds for an off-shell D meson. As it will be seen in the next subsection, the contribution of higher resonances and 
continuum in Eq. (U) will be transferred to the OPE side. 



fofpfD'^mD-mp 
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B. The OPE side 



In the OPE or theoretical side each meson interpolating field appearing in Eqs. ([TJ and ([2]) is written in terms of 
the quark field operators in the following form: 

jP^ix) = d(x)^^u(x)- (10) 



3 (y) = Myhsciy) (ii) 

and 

jr(0) = J(0)7.c(0) (12) 

where u, d and c are the up, down and charm quark field respectively. Each one of these currents has the same 
quantum numbers of the associated meson. The correlators ([T]) and ([2]) receive contributions from all terms in the 
OPE. The first (and dominant) of these contributions comes from the perturbative term and it is represented in Fig.[T] 
Here we will consider the perturbative diagram and the quark condensate. We can write F^^ in terms of the invariant 



D° 




FIG. 1: Perturbative diagrams for the D ofF-shell (left) and p ofF-shell (right) correlators. 



amplitude: 



(13) 



where the meson M(= _D , p) is off-shell. We can write a double dispersion relation for A, over the virtualities p^ and 
p' holding fixed: 



(s, u, t) 



%u-p") 



A 



(M) 



(14) 



where t = and p*^*^^ (s, u, t) is the double discontinuity of the amplitude A'*^^ {p^ , p'"^ , q^) when the meson M{= D , p) 
is off-shell. The perturbative contribution to the double discontinuity in (|14l) for an off-shell D meson is given by: 



{u,s,t) 



71 



M(2m^ — s — t + u) 
A 



(15) 



with A=: {u + s ~ t)'^ — Aus. The integration limits in the integrals in (fT4|) are: 

TO^(s + t — mj) — st 



and 

< s < sq 

Evaluating the perturbative contribution for the double discontinuity for an off-shell p meson we find: 



3toc t 



p^pHu,s,t) = [u + s-t 



2mi 



(16) 
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and the corresponding integration limits in (I14[) are: 

m^(s — t 



< u < Uo 



and 



< s < So 



As usual, we have already transferred the continuum contribution from the hadronic side to the QCD side, through 
the introduction of the continuum thresholds sq and uq . In doing so we made the assumption that at very large 
values of s and u the double discontinuity appearing in the phenomenological side coincides with that of the OPE 
side. This assumption is often called quark-hadron duality. 

In order to improve the matching between the two sides of the sum rules we perform a double Borel transformation 



[16| in the variables = -p^ -)> and P'^ = -p'^ M'- 

Incidentally, this double Borel transform will kill the contribution of the quark condensate a|J!^^ leaving only A^^^* 
which is represented in Fig. [2] and is given by: 



on the invariant amplitude Kqpe and also on Aphen- 



A 



(17) 



where {qq) is the light quark condensate. 



DO 




D*+ im) p+ 

FIG. 2; Contribution of the qq condensate to the D off-shell correlator. 



The sum rule 



After performing the Borel transformation [l^ on both invariant amplitudes A'qp^ and A^^^'^ we identify p4)) 



with ([5]) and then with In doing so we obtain two equations (the sum rules) for the form factors QjjJjjpiQ'^) and 
9D*DpiQ^) respectively. We get the following sum rules: 



C 



[q^ — m?) 



and 



c 



(g2 - ml) 

2 r 

where C = mpfpniD* fo* 



1 

4^ 



1 

4^ 



ds I du p^''\u, s^t) e M^e m' 



ds / dup'^)(u,s,i)e"^e"IF7 + (gg) e"™'/*^' 



(18) 



(19) 



III. RESULTS AND DISCUSSION 

Table |T] shows the values of the parameters used in the present calculation. We used the experimental value for /p 
[l7| , and took /d and fo* from Ref. @. The continuum thresholds are given by sq — (^Tii + As)^ and uq — (toq-I- A„)2, 
where rrii and TOq are the masses of the incoming and outgoing meson respectively. 
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m<:(GeV) 


mD,(GeV) 


mD(GeV) 


mp(GeV) 


/i5*(GeV) 


/D(GeV) 


/p(GeV) 


{qq}{GeYf 


1.27 


2.01 


1.86 


0.775 


0.240 


0.170 


0.161 


(-0.23)^ 



TABLE I: Parameters used in the calculation. 
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FIG. 3: g'oJopiQ^ = 1-0 GeT/^) as a function of the Borel mass 



In this work we use the following relations between the Borel masses M2 and M'2: = for a D off-shell 

and = ^IL^ for a p off-shell. 

Using As — 0.5 GeV and A„ = 0.70 GeV for the continuum thresholds and fixing — 1 GeV^, we found a sum 
rule for g'j^Jjjp as a function of which is very stable with respect to in the interval 20 < A/ ^ < 50 GeV^ . This 
can be seen in Fig. [31 In what follows we choose the value = 30 GeV^ as a reference. In Fig. Uwe show the 
dependence of the form factor g^^ljjp- Here the threshold parameters were taken to be = A„ = 0.5 GeV. Also in 
this case we find a good stability for a wide range of M'^ values. We have chosen the Borel mass to be AP = 3 GeV^. 
Having determined Af^, we calculated the dependence of the form factors. We present the results in Fig.[5l where 
the squares correspond to the gD*^DpiQ^) form factor in the interval where the sum rule is valid. The triangles are 

the result of the sum rule for the ff^*£)p(Q^) form factor. 

In the case of an off-shell D meson, our numerical results can be fitted by the following monopolar parametrization 
(shown by the solid line in Fig. [5]): 

jn) , 234.38 

where the function g'jji^DpiQ^) has the units of GeV~^, as we could anticipate from (jS]). Following our previous works 
0-0, 0, we define the coupling constant as the value of the form factor at = —m\,j, where iriM is the mass of 
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■ Perturvative Contribuition 
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FIG. 4: {('"'(Q^ = 1 GeV^) as a function of the Borel mass. 



the off-shell meson. Therefore, using = ^"Id in Eq ([^0]) . the resulting coupling constant is g^.jjp — 5.76 GeV~^. 

For an off-shell p meson our sum rule results can be fitted by an exponential parametrization, which is represented 
by the dot dashed line in Fig. [SJ 

9^^\u,{Q') = 5.12^-^3^4.33 (21) 

Using Q2 ^ _to2 ^ ^(p)^^ ^ 5 89 QeV-^ 

Looking at Fig. [S] we can observe that the D off-shell form fator is much harder than the p off-shell one. This 
agrees with the conclusions found in most of our previous works: the heavier is the off-shell meson, the harder is its 
form factor. Every extrapolation introduces some ambiguity in the final results, since we have the freedom to fit a 
set of points with different parametrizations. In our case this freedom is strongly reduced because we require that both 
parametrizations lead to the same coupling constant. In Fig. [5] this requirement forces the two endpoints of (j20p and 
(PT|) . which are taken at the squared masses of the corresponding particles, to coincide, i.e., to have the same height 
in the figure. 

In order to study the dependence of our results with the continuum threshold, we vary As,„ between 0.4 GeV < 
As,u < 0.6 GeV in the sum rule (HH) and 0.4 GeV < < 0.6 GeV and 0.65 GeV < A„ < 0.75 GeV in the sum rule 
pop . This variation produces new sets of curves which are shown in Fig. [S] and gives us an uncertainty range in the 
resultig coupling constants. They are g^^^jp = 5.70 ± 0.06 GeV^^ and g^^ljjp = 5.90 ± 0.08 GeV""^. We can see that 
the two cases considered here, off-shell D or p, give compatible results for the coupling constant. Considering the 
uncertainties in the continuum thresholds and taking the average between the obtained values we have: 

go-Dp = (5.8I ± 0.17) GeV"^ (22) 

Our results were obtained for a concrete choice of currents, Eqs. (fTO|) . (fTTj) and (|12l) . which represent charged states. 
Consequently the obtained couplings are for charged states and from them we can get the generic coupling appearing 
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FIG. 5: g'^jtop (squares) and g^j^lj^p (triangles) QCDSR form factors as a function of Q^. The dot-dashed and solid lines 
correspond to the exponential and monopole parametrizations respectively. 



in the Lagrangian ([3]) through the relation: 



go* Dp — 



V2 



V2 



(23) 



Therefore the value of the coupling constant is: 

go-Dp = (5.8I ± 0.17) /V2 = (4.13 ± O.OS) GeV"^ 

In Table II we compare this value with others obtained in previous works. For us the comparison between our 
results and those found in Ref. [1] and Ref. [l^ is especially meaningful, since both approaches use QCD sum rules, 
although in a different implementation. As it can be seen in Table II, these two works arrive at somewhat different 
values of the coupling constant, which are, within the errors, compatible with each other. We use the standard SVZ 
sum rules and the authors of [1, [3 work with QCD Light Cone Sum Rules (LCSR). We use the three-point function, 
whereas they use the two-point function with the p as an external field. The advantage of using the three-point 
function is that it allows us to treat the p meson as an off-shell particle and compute not only the coupling constant 
but also the form factor. Our results have non-perturbative corrections coming from condensates whereas in [l^ 
the authors perform a twist expansion. In view of these differences it is reassuring to see that we obtain values of 
gD'Dp which are compatible with each other. 



This work 


LCSR fl8] 


LCSR fS] 


VDM fl5] 


SU{A) ^ 


4.11 ±0.08 


4.17 ± 1.04 


3.56 ±0.6 


2.82 ±0.1 


3.28 ±0.1 



TABLE II: go-Dp in GeV ^ obtained in previous works. 

In Ref. [T^ the authors made an estimate of the D*Dp coupling constant applying the Vector Dominance Model 
(VDM) to the radiative decay D* — >■ and using experimental information. The obtained value is somewhat smaller 
than the others. We should take this estimate with caution, since it has been known since long ago [l^ that the 
application of VDM to the charm sector is not always reliable. 
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FIG. 6: Dependence of the form factor on the continuum thresholds. The sohd curve corresponds to As,u 
dashed one to As,„ = 0.6 GeV and the dotted curve to As,„ = 0.4 GeV. 



0.5 GeV, the 



Another way to estimate unknown charm coupling constants is to connect them with known couplings through 
SU{A) relations. In the present case, we could use the relation: 



9D'Dp 



go'D' 



6.6 ±0.3 



(3.28±0.15) GeV" 



(24) 



mo* 2.01 

This number is smaller the QCDSR results. In our previous works [1, [l^ we found that, in QCDSR, the SU{4) 
relation gj/^D*D* — 9j/^dd is satisfied. However, from and [3] we observe that other SU{A) relations, such as 

gpD'D' — 9pDD and QpO'D* — ■^gj/4iD*D* are violated at the level of 50 %. This is not surprising since the mass 
difference starts to play an important role when we go from the heavier vector mesons to p. 

In conclusion, we have calculated the form factors of the D*Dp vertex and also the coupling constant. We have used 
QCD sum rules to explore the properties of the three-point Green function of this vertex. The form factors gij^'upiQ^) 
and g^iy>Dp{Q^) (Ell) were obtained for the first time and, as mentioned in the introduction, they can be used in 



several phenomenological applications. The coupling constant extracted from the form factors is go-Dp 
GeV^^ and it is in agreement with other QCDSR estimates. 



4.1 ±0.1 
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